Current-driven orbital order-disorder transition in LaMnOs 

Parthasarathi Mondajj and Diptcn Bhattacharyet] 
Nanostructured Materials Division, Central Glass and Ceramic Research Institute. CSIR. Kolkata 700032. India 



o 

(N 

< 






O 

o 



> 

ON 

rn 

in 

o 
o 



X 



p. Mandal 
Experimental Condensed Matter Physics, Saha Institute of Nuclear Physics, Kolkata 700064, India 

(Dated: August 10, 2011) 

We report a significant influence of electric current on the orbital order-disorder transition in 
LaMnOs. Tlie transition temperature Too, thermal hysteresis in the resistivity (p) versus tem- 
perature(r) plot around Too, and latent heat L associated with the transition decrease with the 
increase in current density. Eventually, at a critical current density, L reaches zero. The transition 
zone, on the other hand, broadens with the increase in current density. The states at ordered, 
disordered, and transition zone are all found to be stable within the time window from ~10^'^ to 



PACS numbers: 71.70.Ej, 64.60.Cn, 71.30.+h, 71.27.-l-a 



I. INTRODUCTION 



The long-range orbital order in LaMnOa develops with 
ordering of active Mn 3d3j.2_r2 and 3d3y2_^2 orbitals, al- 
ternately at each Mn site, within a Mn-0 plane and the 
stacking of this order along c-axis (d-type order) ji It un- 
dergoes a reversible order-disorder transition at a charac- 
teristic transition temperature Too^ The orbital order 
superstructure originates from cooperative fluctuations 
of the doubly degenerate Mn Sde^ orbitals interacting 
via Kugel-Khomskii superexchange. This is further aided 
by cooperative Jahn- Teller distortion of the Mn'^"'"06 
octahedral^"— The structurally forbidden orbital Bragg 
peaks could be clearly observed, with expected azimuthal 
angle dependence of peak intensity, in resonant x-ray 
scattering experiment, thus offering decisive evidence for 
formation of orbital order superstructure in LaMnOai^ 
The orbital order, of course, is not a continuum but 
contains domains due to interaction with intrinsic lat- 
tice strain and/or defects i^ This orbital domain struc- 
ture could also be observed in spatially resolved coherent 
x-ray scattering experiment li It is both technologically 
as well as fundamentally important to explore whether 
or not such an orbital ordered structure undergoes an 
order-disorder transition upon electric, magnetic or opti- 
cal stimulations. It has been shown in the past that the 
long-range charge order in the doped systems melts down 
under electric, magnetic, and optical stimulations yield- 
ing sharp rise in magnetization and/or conductivity to- 
gether with change in the crystallographic structures 4"— 
The orbital stripes too were shown to undergo rota- 
tion under electric field in charge/orbital ordered layered 
manganiteSi ^^'^^ While in LaMnOa, the Jahn- Teller dis- 
tortion and orbital order was found to quench completely 
under ~20 GPa mechanical pressurepi in LaVOs, the 
high energy laser pulse could melt the orbital orderi^^ 
In spite of such rich background, there is, as yet, no in- 
formation about whether or not in undoped LaMnOa the 
orbital order-disorder transition can be driven by an elec- 
tric field or photoirradiation. More importantly, in none 



of the past experiments on external-stimulation-driven 
phase transition, has an attempt been made to track the 
evolution of the order of transition, i.e., how the latent 
heat of the transition, if any, varies with the increase 
in electrical, magnetic, optical, and mechanical energy, 
directly by calorimetry. 

In this paper, we report the observation of significant 
influence of an electric current on orbital order-disorder 
transition in LaMnOs. Using the experimental data, a 
phase diagram on the current density (J)-temperature 
(T) plane is constructed. We find that the latent heat 
(L) associated with the transition becomes zero at a crit- 
ical current density Jc ^ 50 A/cm^. It is also observed 
that the transition width broadens continuously as J is 
increased. The zones in the J-T plane - at well below 
transition, within the transition region, and at well above 
the transition - are characterized by probing their resis- 
tivity relaxation behaviors together with the resistivity 
(p) versus temperature ( T) patterns. 



II. EXPERIMENTS 

The experiments are carried out on high quality sin- 
gle crystals of LaMnOs of dimensions 5 x 3 x (0.5- 
1.5) mm'^J^ The gold electrodes and wires are used in 
standard linear four-probe configuration for the measure- 
ments of p — T and p versus time (at a given tempera- 
ture) patterns under different bias currents. During the 
measurement of resistivity, the temperature sensor is at- 
tached directly onto the surface of the sample in order to 
record the data as a function of actual sample temper- 
ature. We also record the differential thermal analysis 
(DTA) data simultaneously with p — T measurements 
under varying bias current in order to estimate the la- 
tent heat of the transition, from the DTA peak around 
the transition point, as a function of applied current den- 
sity. We report all the data as a function of the applied 
current density as the measurements are done by passing 
different bias current directly through the sample. 
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FIG. 1. (color online) (a) The resistivity (p) versus temperature (T) plot under different bias current at above room temperature 
showing the variation in Too with the bias current density; inset: the region around Too is blown up; (b) The dln{p/T)/d{l/T) 
versus 1/T plots for different bias current density obtained from the resistivity (p) versus temperature (T) data; the temperature 
corresponding to the peak is the Too; inset shows how the Too shifts with the increase in current density in low current density 
regime; (c) ln{p/T) versus 1/T plots; the onset of the transition (T*) is marked by arrow; (d) comparison of the data obtained 
from continuous dc and pulsed current measurement is shown for two representative cases. In all the plots the current density 
increases as one moves from top to bottom. 



III. RESULTS AND DISCUSSION 

In Fig. 1(a) wc show the p -T data measured un- 
der different bias current. It is quite evident from the 
inset of Fig. 1(a) that the transition temperature Tqo 
decreases under increasing current (or bias voltage). In 
Fig. 1(b) we plot the d\n{R/T)/d{l/T) versus 1/T pat- 
terns. The peak in this plot marks the Tqo] the height 
of the peak decreases while the width increases as the 
bias current increases. Figure 1(c) shows the deviation 
of the R-T pattern from adiabatic small polaron hopping 
model R = RqT". exp(A/fcBT) {a = 1) beyond a certain 
temperature T* (marked in the figure by arrow) below 



Too- T*, therefore, marks the onset of the transition. 
An error of ±0.5% is, of course, estimated to be involved 
in identifying the T* . Like Too, T* also decreases pro- 
gressively with the increase in bias current density. The 
zone confined within T* and Too rnarks the transition 
zone where both orbital ordered and disordered phases 
are expected to coexist and the transport of charge carri- 
ers does not follow any model applicable to motion with 
long range coherence. An important issue here is the 
Joule heating of the samples which renders the determi- 
nation of the current driven intrinsic effects difficult. In 
order to quantify the impact of Joule heating under en- 
hanced bias current (or electric field), we measure the rise 



in temperature due to heating by attaching a tempera- 
ture sensor directly onto the sample surface. Such an ar- 
rangement has earlier been used by others for measuring 
the actual sample temperature governed both by heat- 
ing from the bath and Joule heatingi^ We found, as ex- 
pected, that the rise in temperature, as a result of current 
flow through the sample and sample-current lead junc- 
tions, enhances with the enhancement of current density 
J: from nearly negligible to as high as 100 K for a current 
range 50 /xA - lA. We start the current flow through the 
sample at room temperature and allow the temperature 
to rise by a certain extent within that range. Once the 
temperature stabilizes at a particular point above room 
temperature, we start the measurement. At that point 
no difference between the furnace (i.e., bath) and actual 
sample temperature could be noticed. The difference re- 
duces to zero as the sample resistance and, therefore, the 
Joule heating drops down drastically. The actual sample 
temperature has all along been monitored while recording 
the p — T data. Similar dc current driven measurements 
on charge ordered compound has earlier been carried out 
by others and it was found that the impact of Joule heat- 
ing is minimumJ^ In order to quantify the influence of 
Joule heating even further, we have compared the dc re- 
sistivity pattern around Tqo with the pattern obtained 
under pulsed current (pulse width ^500 /is). Both the 
data set are found to be nearly identical in the high tem- 
perature regime (Fig. Id) establishing negligible role of 
Joule heating in that zone. Moreover, reproducibility 
of the features of transition, such as, broadening of the 
peak in the DTA pattern, decrease in the hysteresis and 
jump in p — r around the transition, close matching of 
data between heating and cooling cycles etc in crystals 
of different thickness reveals that these are intrinsic field- 
driven effects. The experiments have been repeated on 
crystals of different thickness and all the features of tran- 
sition were found to be reproducible. The Joule heating 
cannot give rise to these features, reproducibly, around 
Too- 

Figure 2 shows the representative raw DTA thermo- 
grams recorded under different bias current. With the 
increase in current, the peak broadens while the area un- 
der the peak decreases. Since the Joule heating near the 
transition zone is negligible, with the increase in J, the 
baseline of the DTA trace does not change at alli^ There- 
fore, no compensation was necessary for detecting the 
peak and its variation under enhanced current density. 
The latent heat has been calculated by subtracting the 
background by an appropriate technique and identifying 
the peak area properly. The errors in such estimation are 
calculated to be ±0.5% in the case of transition temper- 
atures and ±10% in the case of latent heat. In Fig. 3(a) 
we show the variation in the latent heat {L) of transition 
with J while in Fig. 3(b) the variation of the transition 
temperatures is shown. There is a slight history effect 
as the transition temperatures and peak area differ a bit 
between heating and cooling cycles. This is not because 
of any intrinsic effect (e.g., due to slower phase transition 



dynamics or metastability) as discussed later, but could 
be due to slight impurity in the inert atmosphere (flowing 
nitrogen) maintained during the experiment. In fact, ap- 
plication of different heating/cooling rates did not result 
in any significant shift in the transition temperatures. 
The L is found to decrease gradually with the increase 
in J and reach zero at Jc ~50 A/ cm? [Fig. 3(a)]. This 
pattern of variation of L with J is consistent with the 
variation of the extent of hysteresis (AT) in the p — T 
plot around Too between the heating and cooling runs; 
AT too decreases with the increase in J and vanishes at 
Jc [Fig. 3(a), inset]. The transition zone (marked by 
the onset, peak, and end temperatures) observed in Fig. 
3(b) also broadens with J. It is important to mention 
here that the orbital order-disorder transition even in a 
very high quality single crystal of LaMnOs under nearly 
zero J is actually a broadened first order transition.— No 
thermodynamic evidence for strictly first order transi- 
tion has so far been reported. In comparison, compelling 
thermodynamic evidence for the first order transition has 
been gathered using local magnetization measurement in 
the case of vortex lattice melting in high-To supercon- 
ductor. A step-like rise in magnetization could be noticed 
within a temperature range of ^ 3 mK around the vor- 




720 730 740 750 760 770 780 790 
T(K) 

FIG. 2. The representative raw DTA traces in heating cycle, 
observed under different bias current, showing endothermic 
peaks around the orbital order-disorder transition; the peak 
shifts and broadens while the area decreases as the bias cur- 
rent increases. 



tex lattice melting line<2^ It has also been shown that 
disorder can broaden the first order transitioni^ In the 
present case, of course, the broadening of the transition 
even under very low current density could be because of 
intrinsic inhomogeneity or disorder due to the presence of 
orbital domains. Therefore, we compare the latent heat 
as a function of J only in the sense of noting the relative 
variation. The nature of the transition actually broad- 
ens progressively and finally becomes broader than the 
resolution of the instrument. At that point, the isolation 
of the peak area from the baseline is no longer possible 
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FIG. 3. (color online) (a) The variation of latent heat- 
estimated from the area under the peak observed in DTA 
thermogram near the orbital order-disorder transition point- 
with bias current density (J); inset shows how the extent of 
hysteresis decreases with the increase in bias current density; 
current density increases as one moves from top to bottom 
of the inset plot; (b) variation of the transition temperatures 
(noted from DTA thermograms) with J; obviously the transi- 
tion width increases with J. 



and the estimated L reaches zero. Using a calorimeter of 
higher sensitivity or a local calorimeter one could possi- 
bly resolve the peak. 
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FIG. 4. (color online) The phase diagram of electric current 
driven orbital order-disorder transition on the J-T plane: the 
transition width increases with the increase in J. The peak 
width of the DTA data has been superimposed on the phase 
diagram. It shows that the DTA peak width is quite smaller 
in comparison to the transition width identified from the elec- 
trical resistivity data. For clarity, the error bars have been 
omitted here. The points at which the data of relaxation of 
resistivity have been presented in Fig. 5 are shown as 1,2,3, 
and 4. 

Using the data of transition temperatures as a func- 
tion of bias current, obtained both from electrical and 
calorimetric measurements, we construct a phase dia- 
gram on the J — T plane (Fig. 4). The phase diagram 
indicates three regions: (i) ordered, (ii) disordered, and 
(iii) transition. It is clear from the figure that the tran- 
sition zone broadens progressively with current density. 
Interestingly, the peak broadening, observed in DTA data 
under enhanced J, covers a rather narrow portion of the 
T* — Too transition zone identified from the electrical 
resistivity data. It shows that though the electrical mea- 
surement senses the onset of transition at T* , the calori- 
metric measurement senses the onset at a much higher 
temperature closer to the thermodynamic Too- Earlier 
work2^ on evolution of crystallographic structure across 
T* — ToOi on the contrary, reveals that the anomalous 
structural distortion sets in at ^ T* itself. The reason 
behind discrepancy between the onset points identified 
from the crystallographic and resistivity data and those 
recorded from the calorimetric data could be the differ- 
ence in sensitivity of the probes. The calorimeter sen- 
sitivity is <1 /iW whereas the sensitivity of the nano- 
voltmeter used for recording the voltage drop and hence 
resistance of the sample is <10 nV. Therefore, while 
the electrical resistivity and x-ray diffraction experiments 
could sense the nucleation of the orbital disordered phase 
and hence record accurately the onset of transition, the 



calorimeter could record the onset of transition only when 
the orbital disordered phase has grown beyond a critical 
size and hence at a higher temperature. 
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FIG. 5. (color online) A representative plot of time depen- 
dence of normalized electrical resistance at different tempera- 
tures under a bias current density ~27 A/cni^. The temper- 
atures and the bias current are chosen in such a way so that 
different regions of the phase diagram - ordered, disordered, 
and transition - can be accessed. The extent of relaxation 
appears to be negligible in all the cases indicating stability of 
the phases. On the contrary, 2-4% decay in resistance could 
be noticec^ in the case of metastable orbital ordered phase 
in nanoscale (< 20 nm) LaMnOs. 

In order to probe further the characteristics of these 
three regions of the J — T phase diagram, especially, 
whether or not the slight history effect observed in DTA 
data around the transition is due to intrinsic metastabil- 
ity of the phases at the transition zone, we have measured 
the relaxation of resistivity. For crossing the boundaries 
of the transition zone along a constant J line, we selected 
four points (shown in Fig. 4) with different temperatures. 
The relaxation measurements were carried out by raising 
the temperature of the crystal to the desired point and 
then applying the requisite current density. We reached 
the points separately by raising the temperature and field 
from room temperature and zero-field. The resistivity 
data were recorded over a time span of ~50ms to 3600s 
at an interval of ~50ms after stabilizing the sample tem- 
perature at a given point within less than 0.01 K. We 
repeated such measurements along other constant J lines 
too. In Fig. 5, we show the relaxation characteristics 
observed at points 1, 2, 3, and 4. The characteristics 
are representative of all the points at which similar re- 
laxation measurements were carried out. We observe vir- 
tually no time dependence of resistivity corresponding to 
the points 1, 2, 3, and 4 signifying stability of not only or- 
bital ordered and disordered phases but also of the mixed 
phase within the transition zone. Therefore, the states at 
transition zone are not metastable. The slight difference 
observed in the transition temperatures between heating 



and cooling cycles (Fig. 3) cannot be due to metastabil- 
ity and consequent emergence of transient states around 
the transition. The transition is thermodynamic and the 
transition kinetics is not accessible within the laboratory 
time window^ 10^"^— lO^s at all the temperatures. Prob- 
ing of evolution of crystallographic structure^ within 
and around the transition zone shows that O' orthorhom- 
bic structure (c/v2 < a < 6; a, b, c are lattice pa- 
rameters) of the orbital ordered phase evolves into a 
mixed O' and O orthorhombic phase within the tran- 
sition zone and finally into a pure O orthorhombic phase 
(c/v2 ~ a ~ &) above the transition zone. By taking 
into consideration the evolution of the structure within 
the transition zone along with the relaxation data of re- 
sistivity, it is possible to conclude that the mixed 0'-|-0 
phase does not have any temporal fiuctuation within the 
laboratory time window. On the contrary, orbital or- 
dered phase in nanoscale (<20 nm) LaMnOa is found to 
be metastable with 2-4% decay in resistance with time 
and irreversible order-disorder transition within similar 
time scale<^ Emergence of transient metastable states 
with lifetime (~ 10""^ — lO'^s) comparable to the labo- 
ratory time scale has been observed by Kalisky et al^ 
around the vortex solid-solid phase transition as well, in 
a high-Tc superconductor. 

The decrease in Too Si,s well as latent heat and hys- 
teresis associated with the transition possibly result from 
electric current driven depinning of orbital domains. The 
orbital ordered phase even in stoichiometric LaMnOs is 
not a continuum under zero electric current. It con- 
tains domains due to interaction with intrinsic elec- 
tronic/lattice defects, strain etd^ii Like in the case of 
charge density waves in solids?^ these intrinsic defects 
act as pinning centers for orbital domains as well. In fact, 
orbital order-disorder transition takes place via depin- 
ning of orbital domains.— In the orbital disordered state, 
the short-range order, with high mobility and hence tem- 
poral fluctuations, prevails.— Application of enhanced 
current density (bias field) leads to electromigration of 
defects^ which, in turn, can give rise to field driven de- 
pinning transition. The depinning transition of charge 
density waves as well as vortex lattice in high-Tc super- 
conductors has been thoroughly studied 1^ It has been 
shown that depending on the concentration of pinning 
centers and applied force, the depinning of the charge 
density waves can either follow a two-stage or a single- 
stage process. If the concentration of defects is strong, 
the domains start sliding plastically under a small force 
which finally gives way to a sharp transition into a co- 
herent collective movement in 3D under higher force. On 
the other hand, in the case of weak disorder or low con- 
centration of defects, the depinning transition becomes a 
continuous process and yields a coherently moving col- 
lective state continuously. The weak disorder model pre- 
dicts that the depinning transition yields an exponential 
variation of the transition energy scale (e.g., transition 
temperature) with the applied forcci^ Interestingly, this 
model is found to be valid in the present case. It has been 



observed that the Too versus current density J pattern 
(Fig. 4) fohows closely the model Too(J)/Too(0) = 
exp[-J/Jo] (dashed line, Fig. 4), except at a very low 
J (Jo is a constant here and Too is the transition tem- 
perature under zero field). Validity of this model in the 
present case provides an indirect support for the conjec- 
ture of current density (bias field) driven dcpinning of 
orbital domains. Because of variation in the depth of the 
potential well of defects, the screening of direct force for 
electromigration will vary^l which, in turn, is expected 
to give rise to inhomogeneous dcpinning. This inhomo- 
geneous depinning is possibly the origin of broadening of 
the orbital order-disorder transition zone and drop in the 
latent heat. 

Can there be any other origin for the influence of 
electric field on orbital order in LaMnOa? The d-typc 
orbital order in LaMnOa cannot give rise to a elec- 
tric dipole moment. It only produces a higher order 
quadrupole moment^ which cannot couple linearly with 
the applied electric field. Even the domain boundaries 
of orbital order cannot be intrinsically charged. There- 
fore, unlike charge order, microscopically, long-range or- 
bital order in LaMnOs should not be influenced by elec- 
tric field. Whether unleashing of charge carriers via 
Mn^~^ — > Mn?^ + Mn^^ disproportionation reaction un- 
der field could then be the origin of such an effect? The 
generation of mobile charge carriers under field would 
have given rise to even more dramatic effect. There- 



fore, it seems, apparently, that these effects are not really 
playing any significant role here. Study of orbital do- 
main structure under external electric field using spatio- 
temporally resolved resonant x-ray scattering data can 
offer direct proof of the orbital domain depinning transi- 
tion under field. This is beyond the scope of this paper. 



IV. SUMMARY 

In summary, we observe significant influence of electric 
current density (bias field) on orbital order-disorder tran- 
sition in a single crystal of pure LaMnOa . The transition 
temperature, hysteresis, and the latent heat of transition 
decrease monotonically with the increase in the current 
density; the transition width, on the contrary, increases. 
Finally, at a critical current density, the latent heat be- 
comes zero. The states at the ordered, disordered, and 
transition zones are all found to be thermodynamically 
stable within the laboratory time scale. This current 
density driven orbital order-disorder transition possibly 
originates from the field driven depinning of orbital do- 
mains via electromigration of intrinsic defects. 
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